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Abstract. Canonical Discriminant (CDA), ANOVA and AN0S1M analyses were calculated for all recently known dis¬ 
tribution areas of Iberolacerta cyreni including several small and unstudied peripheral populations. The only differenti¬ 
ated sample is Guadarrama (the nominate subspecies), with very limited overlap in the CDA (correct classification > 70%) 
and different from nearly all the other samples in ANOSIM. Guadarrama is a recently differentiated but well diagnos- 
able (morpho)subspecies (with lower values of dorsalia, ventralia and greater values of circumanalia). Despite the mtD- 
NA differences of the Bejar specimens, their morphology is largely equivalent to that of I. cyreni castilianci (Gredos), 
but clearly differ in their female body elongation (near 1 cm) with shorter limbs, a possible strategy to increase clutch 
size. Populations from the Sierras de Avila (Villafranca, Serrota and Paramera) are very similar among them. Villafran¬ 
ca (in males) together with Bejar (in females) are the most connected samples in MST, and the root of the species dif¬ 
ferentiation from a morphological point of view, once discarded geographical and climatic influence on morphology. All 
populations except Guadarrama shall be considered as/, c. castiliana by their morphological identity with Gredos. These 
morphological similarities probably are the reflect of extensive gene flow among them, responsible of maintaining their 
morphology largely equivalent. 

Key words. Lacertidae, Iberolacerta cyreni , Intraspecific variability, subspecies, Geographical variation, Iberian Penin¬ 
sula. 

Resumen. Se ha calculado un Analisis Discriminante Canonico, ANOVA y ANOSIM con toda el area de distribucion de 
Iberolacerta cyreni , incluyendo varias poblaciones perifericas no estudiadas hasta la fecha. La unica muestra diferencia- 
da es Guadarrama (la subespecie nominal), con niuy poco solapamiento en el CDA (clasificacion correcta > 70%) que 
difiere de practicamente todas las demas muestras en el ANOSIM. Guadarrama es una poblacion recientemente diferen- 
ciada, pero bien diagnosticable como (morfo)subespecie (valores bajos de dorsalia y ventralia, y altos de circumanalia). 
A pesar de las diferencias mitocondriales de Bejar, su morfologia es ampliamente asimilable a /. c. castiliana (Gredos), 
siendo destacable el relativo elongamiento coiporal de las hembras (casi 1 cm) con miembros proporcionalmente cortos, 
una posible estrategia para incrementar el tamano de puesta. Las poblaciones de las Sierras de Avila (Villafranca, Serro¬ 
ta y Paramera) son muy similares entre si. Villafranca junto con Bejar (en machos y hembras respectivamente) estan mor- 
fologicamente en la raiz de la diferenciacion de la especie (MST), una vez descartada cualquier influencia climatica o de 
distancia geografica. Excepto Guadarrama, todas deben considerarsc como I. c. castiliana por su identidad morfologica 
con Gredos, lo que reflcja la probable prcscncia de un flujo genetico extensivo y recientc entre ellas. 


INTRODUCTION 

The Spanish Sistema Central consists of a series of Sier¬ 
ras, more or less aligned in a ENE-WSW direction, which 
separate the Duero (to the North) and Tajo (to the South) 
river drainages, or what is the same, the Old and New 
Castile plateauxes. It runs from the Portuguese Serra da 
Estrela (inhabited by a relict population of Iberolacerta 
monticolci [Boulenger, 1905]), across the Spanish Sierra 
de Gata (apparently too low and dry for Iberolacerta), the 
Sierra de Francia witli Iberolacerta martinezricai (Arribas, 
1996) and the main pail of the Spanish portion of its range 


(over 240 km in length), which is inhabited in several 
points by populations of the Carpetane Lizard ( Iberolac¬ 
erta cyreni [Muller & Hellmich, 1937]). 

The Carpetane Lizard is widespread through the main 
parts of the Sistema Central and is mainly known from 
Sierras de Bejar, Gredos and Guadarrama (Fig. 1). It was 
raised to species level (Arribas 1996) based on allozymes 
(Mayer & Arribas 1996), karyology (heterochromatiniza- 
tion of sex-chromosome and localization of the NORs; 
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Odierna et al. 1996) and adult and hatchlings pattern and 
coloration. Two subspecies were defined, the nominal I. 
c. cyreni from Guadarrama (type locality: Puerto de 
Navacerrada), and /. cyreni castiliana from Gredos (type 
locality: Circo de Gredos, Avila) to which frequently are 
assimilated Bejar specimens. This latter subspecies differs 
from the nominate one by a reduced dark pattern, more 
dorsalia, ventralia, slightly larger diameter of the masse¬ 
teric and hindlimb length, and lower circumanalia (Arribas 
1996). 

The degree of genetic differentiation between I. c. cyreni 
(Guadarrama) and I. c. castiliana (Gredos) was analyzed 
by Mayer & Arribas (2003), who found a mtDNA se¬ 
quence divergence of 0.6 % in the 12s rRNA (12s) and 
16s rRNA (16s) mitochondrial genes, which corresponds 
to approximately 0.6 MY BP. Carranza et al. (2004) sug¬ 
gested that both subspecies diverged approximately 
0.8±0.2 MY BP, an estimation mainly based on the Cy¬ 
tochrome b (Cyt b) mitochondrial coding gene (the 12s 
and the nuclear gene C-mos were uninformative at this lev¬ 
el), a divergence time almost identical to the one calcu¬ 
lated by Crochet ct al. (2004) using also the Cytb gene [1.6 
% genetic divergence, which roughly corresponds to 0.6 
to 1, with a mean of 0.8 MY BP). These two values were 
very similar to the above-mentioned ciphers. The inferred 
divergence time increased up to 1.2 (Cyt b) or 1.6 MY BP 
(Cyt b+12s) when different terminal taxa evolutionary 
models and phylogenetic methods were used (Arribas et 
al. 2006; Arnold et al. 2007; respectively). 

On the other hand, specimens from Sierra de Bejar 
branched at the base of the /. cyreni clade in some mtD¬ 
NA analyses (Carranza et al. 2004; Arribas & Carranza 
2004). It was suggested that the split of this populations 
occurred approximately 1,7±0.3 MY BP (Carranza et al. 
2004). However, in analyses using the same mtDNA re¬ 
gions but different taxa and other evolutionary models and 
phylogenetic methods than above (Arribas et al. 2006; 
Arnold et al. 2007), the specimens from Bejar fonned a 
trichotomy with /. c. castiliana from Sierra de Gredos and 
I. c. cyreni from Navacerrada. 

Apart from the uninformative C-mos nuclear gene frag¬ 
ment analyzed by Carranza et al. (2004) there is only one 
other infonnation about differences at the nuclear level, 
the analysis of allozyme data by Almeida et al. (2002), 
which showed a Nei’s distance of 0.002 between speci¬ 
mens from Gredos and Guadarrama. 

From West to East, the distribution of I. cyreni can be di¬ 
vided into two axes connected by low mountain valleys 
(see appendix II), but not clearly interrupted by clear cut 
barriers as river valleys. One axis runs across Sierra de 
Bejar (summit in La Ceja, 2.425 m) and Gredos (Alman- 


zor, 2.592 m), whereas the other axis is constituted by the 
Sierras de Villafranca (Moros, 2.065 m), La Serrota (Ser- 
rota, 2294 m), La Paramcra (Zapatero, 2.160 m) and 
slightly separated by lower areas, Guadarrama (Penalara, 
2.430 m). The two axes greatly overlap longitudinally 
leaving the Villafranca, Serrota and Parameras just to the 
North of the Sierra de Gredos, but at their orographic shad¬ 
ow for rains, and climatically more continentalised. This 
explains the botanical similarities between the Paramera- 
Serrota-Villafranca axis and the Sierra de Guadarrama 
(Luceno and Vargas 1991). 

In Guadarrama (where Podarcis muralis also exists), I. 
cyreni occurs only at the highest areas, from 1.760 m 
(Puerto de Cotos, Puerto de Navacerrada) up to the peaks 
(2.340 m in Penalara). In Gredos it lives almost from 1.700 
to 2.500 m. It was seen in 17. Vll. 1986 (own data) in 
Puerto del Pico (at 1.352 m close to one of the fountains 
of the pass) but recent research in this area has been to¬ 
tally unfruitful. It is possible that these lower stations 
favoured by accelerated cold winds in the mountain pass¬ 
es (Venturi effect) had disappeared by climatic or best, by 
habitat degradation due to human over-frequentation dur¬ 
ing the last 20 years. In Bejar it has been found between 
1.837 m (own data) and 2.443 m (see Lizana et al. 1988, 
1992, 1993; and Martin 2005 for general data; own data 
corrections for the confirmed lower limits). 

Apart from the better known Sierras, the presence of the 
Carpetane Lizard in the small parallel mountain ranges 
called “Sierras de Avila” or “Parameras” (composed by 
three Sierras: Villafranca, La Serrota, and La Paramera) 
was first discovered by the mountaineering group “Valle 
de Ambles” (Lizana et al. 1993), but no specimens have 
been studied so far. All aspects of morphology, status and 
relationships of these small and isolated populations from 
the Sierras de Avila are totally unknown. In these Sierras 
the species is extremely localized, especially in La Ser¬ 
rota and Paramera. In Sierra de Villafranca, the area with 
a relatively more extended suitable area, 1 have found it 
from 1.850 m probably up to the highest areas (Pico Mo¬ 
ros, 2.065 m). In La Serrota it is extremely rare and lo¬ 
calized, also cornered in the highest parts, from 2.284 m 
(perhaps 1.935 m where excrements, possibly of this 
species, were seen; pers. obs.) to the very summit (Pico 
Serrota, 2.294 m); and in La Paramera from 1.700 m in 
the northern slopes to the summit (Pico Zapatero, 2.160 
m) (own data). 

After a three-year prospection of these parallel ranges, I 
gathered data from these localized and barely known pop¬ 
ulations in order to check the relationships of all the Car¬ 
petane Lizards throughout its range. My aim is: a) to re¬ 
assess differences between I. c. cyreni and I. c. castiliana 
in the light of the existence of other small and isolated pop- 
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illations: b) to ascertain the taxonomic status of the Bejar 
populations and to check if these represent a further sub¬ 
species; c) to study both the relationships among the sam¬ 
ples from the Sierras dc Avila (=Parameras) massifs, as 
well as their similitude and differences with their neigh¬ 
bouring and well known populations from Gredos, or the 
more distant populations from Guadarrama and Bejar; and 
d) as the type series of I. cyreni was destroyed during the 
Second World War (SWW), to choose a Neotype for the 
species (see appendix I) in order to fix unequivocally the 
type locality (although apparently all lost, there were al¬ 
so specimens from Gredos in the original type series). 

MATERIAL & METHODS 

Morphology 

A total of 106/92 male specimens, and 136/135 female 
specimens of /. cyreni with a complete measurements 
dataset and snout-vent length greater than 45 mm, were 
included in the univariate (ANOVA) and multivariate (dis¬ 
criminant) analyses, respectively. Given that these popu¬ 
lations present sexual dimorphism (Arribas 1996, 1999a; 
Arribas et al. 2006), analyses were carried out for males 
and females separately. All material is from Oscar Arribas 
(OA) database. 

OTUs names, localities and specimens included in the 
morphological multivariate analysis were as follows (Fig. 
1 ): 

GUADARRAMA: Sierra de Guadarrama (Madrid and 
Segovia provinces, Spain). 25 males and 36 females [I. 
cyreni cyreni ]. MIJARES: Puerto de Mijares (Sierra del 



Fig. 1 . Schematic representation of the distribution of Ibero- 
lacerta cyreni in the Spanish Sistema Central. The different lo¬ 
calities (OTUs) cited in the text are represented. 1: Bejar; 2: Gre¬ 
dos; 3: Villafranca; 4: Serrota; 5: Paramera; 6: Mijares; 7: Gua¬ 
darrama. 


Cabezo, Gredos Oriental Massif, Avila province, Spain). 
6 males and 7 females. GREDOS; Circo de Gredos (Gre¬ 
dos Central Massif, Avila province, Spain). 23 males and 
46 females [/. cyreni casliliana], BEJAR: Sierra de Be¬ 
jar (Gredos Occidental Massif, also known as Sierra de 
Candelario, Salamanca province, Spain). 11 males and 28 
females. VILLAFRANCA: Sierra dc Villafranca (Avila 
province, Spain). 20 males and 14 females. SERROTA: 
La Serrota (Avila province, Spain). 3 males and 2 females. 
PARAMERA: Sierra de La Paramera (Avila province, 
Spain). 4 males and 1 female. 

These populations are discontinuous among them (Appen¬ 
dix II) and constitute discrete geographical OTUs (Fig. 1). 

Due to lower sample, it was necessary to cluster the Ser¬ 
rota and Paramera specimens in the male discriminant 
analysis, and these two plus Villafranca (all of them “Sier¬ 
ras de Avila”) in the female one. However, reciprocal dis¬ 
tances between each one of these poorly represented sam¬ 
ples to the best represented ones were carefully checked 
and commented in the results section. As the three popu¬ 
lations from Sierras de Avila seemed to be largely equiv¬ 
alent in the multivariate analyses, a posteriori, all of them 
were treated as a single OTU (S a AVILA) in ANOVA. 

Characters studied 

Biometric characters. Snout-vent length (SVL); Forelimb 
length (FLL); Hindlimb length (HLL); Pileus length (PL); 
Pileus width (PW); Parietal length (PaL); Masseteric scale 
diameter (DM); Tympanic scale diameter (DT); Anal 
width (AW) and Anal length (AL). All linear measure¬ 
ments were made with a digital calliper to the nearest 0.01 
mm. These measurements were transformed to the follow¬ 
ing more informative and not dimensional-depending ra¬ 
tios: FLL/SVL (relative forelimb length; “FLL index”); 
HLL/SVL (relative hindlinb length, “HLL index”); 
PL/PW (pileus shape, “Pileus index”); DM/PaL (relative 
masseteric plate size, “Masseteric index”); DT/PaL (rel¬ 
ative tympanic size, “Tympanic index”); AL/AW (anal 
plate surface, “Anal form index”) and AS/SVL 
(V(AL*AW)* 100/SVL, relative anal plate size with re¬ 
spect to the total length, “Anal size index”) (Arribas 1996, 
2001). The results of the linear measurements and index¬ 
es yielded largely similar results. All ratios were given 
multiplied by 100 to avoid excessive decimal scores. 

Scalation characters. Supraciliar Granula (GrS) for the 
right and left sides; Gularia (GUL); Collaria (COLL); Dor- 
salia (DORS); Ventralia (VENT); Femoralia rigth (FEMr) 
and left (FEM1); 4 th . digit Lamellae (LAM); and Circum- 
analia (CIRCA). 
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Statistical Procedures 

Statistical analyses used in the morphological study in¬ 
cluded both Univariate (ANOVA for SVL, scalation char¬ 
acters and indexes, with post-hoc Tukey-Kramer tests at 
P < 0.05 and P < 0.01 to detect differences among sam¬ 
ples) as well as Multivariate techniques (Canonical Dis¬ 
criminant Analysis, CDA). In this later analysis, each pop¬ 
ulation is represented by a centroid (a hypothetical mid¬ 
dle individual). Minimum-length spanning tree (MST) was 
computed from the Mahalanobis’ distance matrix to de¬ 
tect the nearest neighbours based on their position in the 
multidimensional space. MST representation also avoids 
distortion of UPGMA trees by the reciprocal pairwise dis¬ 
tances recalculation in every step during their construc¬ 
tion. UPGMA frequently clusters samples reflecting sam¬ 
ple sizes than their true relationships. Distances of small 
samples or isolated specimens appear greatly exaggerat¬ 
ed with respect to the well represented ones. As a result 
of that, the small-sized samples appear ever as the most 
external or differentiated in UPGMA derived trees 
(Kunkel et al. 1980; Cherry et al. 1982; Arribas 1997). 
Moreover, the UPGMA trees based in very unevenly sized 
samples also gave very poor Cophenetic Correlation In¬ 
dexes between the tree-derived ultrametric distances ma¬ 
trix and the original Mahalanobis distance matrix and 
therefore we have not used them (Arribas et al. 2006). 

To test the significance of the differences among pre-es¬ 
tablished groups for the Discriminant Analysis (based in 
a geographical origin), we carried out an Analysis of Sim¬ 
ilarity (ANOSIM) (Clark 1988, 1993) that tests if the as¬ 
signed groups are meaningful, this is, more similar with¬ 
in groups than with samples from different groups. The 
method uses the Bray-Curtis measure of similarity to con¬ 
struct clusters of specimens. The null hypothesis is there¬ 
fore that there are no differences between the members of 
the compared groups (they are randomly blended). R-sta- 
tistic scales from +1 to -1. Values closer +1 correspond 
to a perfect case in which all groups were completely dif¬ 
ferent (all specimens of the same group are more similar 
among them than to any specimens of the other groups). 
R = 0 occurs if the high and low similarities are perfect¬ 
ly mixed and bear no relationship to the group, a common 
situation if some of the groups are largely equivalent. A 
value of -1 indicates that the most similar samples are all 
outside of the groups (all groups largely equivalent and 
randomly formed). To check for significance, pseudorepli¬ 
cation tests (1000 randomizations) were run to test if the 
given results can occur by chance. If the value of R is sig¬ 
nificant, there is evidence that the samples within groups 
are more similar than would be expected by random 
chance. Even more important, pairwise tests among 


populations permit to test significance of the differences 
among the concerned groups and to detect which ones are 
really different from the others. 

Mahalanobis’ (squared) distance matrices were compared 
by means of Mantel Test (with 1000 permutations) with 
matrices composed by Euclidean (squared) distances for 
the climatic characteristics of localities: a) Precipitation 
(mm) during the incubation months (July and August, as 
scalation is invariant during lizard’s life); b) Annual pre¬ 
cipitation (mm), c) Temperature (°C) (July and August); 
d) Annual average temperature e) Sun radiation (n 10 
kJ/(m 2 *day*micrometer,)(July and August), and f) Annu¬ 
al Sun radiation. Data were extracted from Ninyerola et 
al. (2005). Also, these Mahalanobis’ distances were com¬ 
pared with (d) the aerial (straight) geographical distances 
among the sampling localities. 

Multivariate (Discriminant and ANOSIM) analyses were 
performed with Community Analysis Package 4.0 (Hen¬ 
derson & Seaby 2007). MST trees and Mantel tests were 
calculated with NTSYS 2.1© (Rohlf2000). Univariate sta¬ 
tistics were processed with NCSS 2001© package (Hintze 
2001 ). 

RESULTS 

Males 

Canonical Discriminant Analysis. The CDA conducted 
with 106 male specimens shows three significant axes that 
explain an 85 % of the total intersample variation. The two 
first axes together explain the main part (70.6 %) and dis¬ 
criminate fairly well the samples, especially the first one, 
the unique with an eigenvalue >1. The first discriminant 
axis has an Eigenvalue of 1.54 (51.2 % of variance ex¬ 
plained; Chi-Sq. with 85 df = 200.71, P < 0.0001) and dis¬ 
tributes the samples with fairly overlap among them, ex¬ 
cept Guadarrama, that has only a small coincidence with 
the other ones (Fig. 2 A). Guadarrama appears in the neg¬ 
ative part of the axis, characterised by the lower values 
for DORS (0.441553) and VENT (0.560765) and greater 
values of CIRCA (-0.572683). Second and third axes 
(eigenvalues < 1) present a considerable overlap among 
the samples and do not discriminate populations. 

This discriminant analysis applied to the samples reached 
a 72.6% of correct classification among the specimens. 
The Guadarrama sample (/. cyreni cyreni ) reaches a 71.9 
% of correct classification in respect to all the other sam¬ 
ples (/. c. castiliana). 
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Discriminant Plot - ib. cyreni MALES 



Discriminant Plot - ib. cyreni FEMALES 



Fig. 2. Canonical Discriminant Analysis (CD A) plots for a) males (above) and b) females (below). Specimens, sample centroids 
and group perimeters are represented. Guadarrama (inverted triangles), Gredos (triangles), Bejar (irregular circles), Mijares (cross), 
Villafranca (asterisk) and Serrota-Paramera (sail). In females, the three last samples are grouped as Sierras de Avila (sail). Sample 
centroids are represented by a square. See text for axis characteristics and results. 


Minimum-length spanning tree (not represented) connect¬ 
ing the centroid (hypothetical middle specimens) of each 
sample is fairly congruent with their geographical posi¬ 
tion, connecting in general neighbouring samples. The 
most “central” (most connected) population is Villafran¬ 
ca that connects with Gredos (at Mahalanobis Distance of 
3.1870), Bejar (3.2199), Paramera (4.8042) and finally, to 
the most isolated one, Guadarrama (6.8412). Two popu¬ 
lations show overestimated distances due to their small 
sample sizes: Mijares (East Gredos) that connects with Be¬ 
jar (7.0951), and Serrota with their neighbouring Param¬ 
era (9.1822). 


Analysis of Similarity (ANOSIM) (Table 3) shows that 
there is a considerable overlap among samples (R-statis- 
tic = 0.122088, P < 0.005; 1000 randomizations) as our 
value (that can range from 1 to -1), although positive, is 
very small. Very significant differences among the (geo¬ 
graphically) assigned groups, appear only among Guadar¬ 
rama and Bejar, Gredos and Villafranca (P < 0.01) but do 
not reach significance with Mijares and Serrota + Para¬ 
mera (both with small samples). The other populations are 
not differentiated among them (P > 0.01). 
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Fig. 3. Iberolacerta cyreni ccistiliana. a) La Covatilla Sky re¬ 
sort (Sierra de Bejar), July 2007, Male ; b) El Travieso (S a de 
Bejar), July 2004, Female ; c) El Calvitero (Sierra de Bejar), Ju¬ 
ly 2004, Female (atypical pattern, with diffumination and coa¬ 
lescence in a unique vertebral line); d) Puerto de Mijares (Gre- 
dos Oriental Massif), July 2006, Female. 


The Analysis of Variance (ANOVA) (Appendix III, Table 
I) indicate that Guadarrama differs from all or nearly all 
the other populations in VENT and CIRCA (with the 
smaller and greater values for these parameters, respec¬ 
tively, in the former population), but also appeared dif¬ 
ferences between Guadarrama and Bejar in Dors (small¬ 
er in the former), and with Gredos in PV (greater in the 
former). An interesting and significant difference appears 
in DORS among Gredos and Bejar samples (clearly 
greater in the later). 

There is no significant correlation among Mahalanobis’ 
distances and any of the geographic and climatic param¬ 
eters analyzed (all Mantel Tests P > 0.05). 

Females 

Canonical Discriminant Analysis: The CDA conducted 
with 136 female specimens shows three significant axes 
that explain a 93.8 % of the total intersample variation. 
The two first axes together explain a large part of the vari¬ 
ance (85.7 %), and especially along the first one, that ac¬ 
counts itself for 62.7 % of the total variation and is the 
unique with an eigenvalue >1 (1.77), discriminating 
Guadarrama specimens from the other neighbour samples 
only with a small overlap (Fig. 2B). The other samples 
show a considerable overlap among them. Guadarrama 
discriminates towards the negative part of the axis, char¬ 
acterised by the lower values of DORS (0.62) and VENT 
(0.66) and greater ones of CIRCA (-0.38). Second and 
third axes (eigenvalues < 1) present a considerable over¬ 
lap among the samples and do not discriminate popula¬ 
tions. 

The discriminant analysis applied to the samples reached 
a 74.26% of correct classification among the specimens. 
Guadarrama sample (I. c. cyreni ) reaches an 87.2 % of cor¬ 
rect classification with respect all the other samples (/. c. 
ccistiliana). 

Minimum-length spanning tree (not represented) connect¬ 
ing centroids is very similar to the male one. The most con¬ 
nected sample is Bejar, which clusters with Villafranca (at 
2.9), Gredos (3.5) and Mijares (7.38, but here exaggerat¬ 
ed by the scarce sample of the later). Guadarrama connects 
with the scarcely represented (and geographically inter¬ 
mediate) Mijares (East Gredos) (at 6.09), and all the Sier¬ 
ras de Avila samples cluster together (Villafranca with Scr- 
rota + Paramera at 8.66). 

Analysis of Similarity (ANOSIM) (Appendix III, Table 
3) shows that there is a considerable overlap among sam¬ 
ples (R-statistic = 0.162588, P < 0.001; 1000 randomiza¬ 
tions), but the results are slightly best than for the male 
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Fig. 4. Iberolacerta cyreni castiliana. a) Pico Zapatero (Si¬ 
erra de la Paramera), July 2005, Male; b) Puerto de Pena Negra 
(Sierra de Villafranca), July 2006, Male; c) Pico Serrota (La Ser- 
rota Massif), July 2005, Female; d) Pico Serrota (La Serrota Mas¬ 
sif), July 2006, Male. 


analysis. Very significant differences (P < 0.01) appear 
among Guadarrama and all the other samples except with 
Mijares. The other populations are not differentiated 
among them (P > 0.01). 

The Analysis of Variance (ANOVA) (Appendix III. Table 
2), as in the male analysis, it shows that Guadarrama is 
the most different one, especially in DORS, VENT and 
CIRCA (the first two characters smaller, and the third one 
greater in the former population). Guadarrama also dif¬ 
fers from Gredos by its lower GUL, from Bejar by its 
greater FLL, HLL, a lower SVL; and from Sierras de Avi¬ 
la by its greater relative anal scale surface. 

Also, significant differences appear in SVL between Be¬ 
jar (clearly the great sized female population) and Gre¬ 
dos, and among this latter (with relative greater FLL and 
HLL) with Bejar and Villafranca. 

As in male analysis, there is no significant correlation 
among Mahalanobis’ distances and the geographic and cli¬ 
matic parameters analyzed (all Mantel Tests P > 0.05). 


DISCUSSION 

From the Discriminant and ANOVA analyses it appears 
that the only differentiated sample is Guadarrama. It ap¬ 
pears with very limited overlap with the other samples in 
the CDA graphs (Figs 2 A and B). Diagnostic characters 
for this population (nominate subspecies: I. c. cyreni) are 
the lower values of DORS (difference more marked in fe¬ 
males), lower values of VENT and greater CIRCA. More¬ 
over, ANOSIM analyses show that Guadarrama is the 
unique OTU that is significantly different from nearly all 
the other samples, except from the close population of Mi¬ 
jares (in both sexes) and Serrota + Paramera (but these ex¬ 
ceptions occur only in the males and probably due to their 
scarce sample size). 

Mijares sample (very small) seems in some aspects ap¬ 
proaching to Guadarrama (specially in DORS and VENT 
values) but globally are clearly closer to I. c. castiliana 
(specially to Bejar in male and female MST). 

Populations West from Guadarrama show a great overlap 
in CDA and lack differences in ANOSIM, being morpho¬ 
logically fairly equivalent and all of them assimilable to 
I. c. castiliana. There are only a few scattered very sig¬ 
nificant differences among them (P < 0.01) in ANOVA, 
as for instance among Gredos and Bejar (this latter has 
greater DORS and a strikingly greater SVL and propor¬ 
tionally shorter limbs that Gredos, but only in female spec¬ 
imens). The reason of the longer SVL in Bejar females 
(from 8 mm to 1 cm greater than in the other populations) 
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which leaves proportionately shorter limbs, can be a con¬ 
sequence of body elongation that in lacertids appears 
linked to a greater clutch size (Brana 1996). This is an in¬ 
teresting question for future study: if Bejar specimens ef¬ 
fectively have greater clutch size that other I. cyreni pop¬ 
ulations. 

Both the MST results (in which Villafranca and Bejar are 
the most connected samples) as well as the presence of 
related species further West (/. martinezricai and I. mon- 
ticola), suggest an origin of the species towards the west¬ 
ern extreme of their current distribution area. From these 
westernmost parts, where it also occurs the higher haplo- 
type diversity (see below), I. cyreni spread towards the 
East. Despite that the Sierras de Avila (Villafranca, Ser- 
rota and Paramera) are slightly more aligned with the 
Guadarrama axis than with Gredos one, we cannot be sure 
from which of these two mountain ranges the former was 
colonized, as MST results in males and females are con¬ 
tradictory. According to the male analysis Guadarrama is 
more related to Villafranca, whereas in the female analy¬ 
ses, it is Mijares (Eastern Gredos) the most related one. 

The results of the mitochondrial analyses of these sam¬ 
ples (Cyt B and 12s) (unpublished, Carranza, pers. com.) 
indicates that the interruption of gene flow is fairly recent, 
as a common haplotype appears in all populations except 
Villafranca and Bejar. All Gredos, Guadarrama and La 
Serrota specimens are identical for these two mitochon¬ 
drial fragments. Independent changes in one nucleotide 
with respect to the common haplotype appear in Villafran¬ 
ca (the unique change is different in two specimens), Para¬ 
mera, Mijares and Bejar specimens, and two changes ac¬ 
cumulate in one Bejar and one Mijares specimen (others 
have only one). 

The current morphological differences of Guadarrama 
specimens seem to be relatively recent, and are possibly 
the result of bottleneck effects during the colonization 
process, or alternatively of strong selective pressures (or 
a combination of both causes). Conversely, the absence 
of marked differences among the other populations 
(more or less with a similar age) could be due to the main¬ 
tenance of a more continuous gene flow among them, re¬ 
sponsible of maintaining their morphology largely equiv¬ 
alent (nuclear genes remain unstudied). The current larg¬ 
er geographical gap in the distribution of I. cyreni occurs 
precisely between Guadarrama and the remaining popu¬ 
lations to the West. 

Despite the presumably short isolation time, as comment¬ 
ed above, a considerable selection pressure or a genetic 
bottleneck in the expanding populations might have pro¬ 
moted and fixed the morphological differences now seen 
in Guadarrama specimens. These factors do not seem to 


be due to isolation-by-distance processes but by histori¬ 
cal vicariant events (cf. Irwin 2002) as there is no rela¬ 
tionship between morphological differentiation and geo¬ 
graphic distances. Also there is no relationship among the 
more obvious climatic parameters (precipitation, temper¬ 
ature and sun radiation) and these differences. 

Concerning the position of Bejar populations, only spec¬ 
imens from the west-facing slopes of the massif have been 
studied, and therefore it is possible that in other parts of 
the massif other haplotypes (the common one with Gre¬ 
dos) could be present. The species was cited from “Puer¬ 
to de Tomavacas, SA” (for Salamanca, sic.!; a mistake as 
this locality is in Avila) (Lizana et al. 1992). This is the 
natural pass between Bejar and Gredos, but I have been 
unable to find it there. This place is a sub-Mediterranean 
environment with Pyrenean Oak open forest inhabiting 
populations of Timon lepidus, Psammodromns algirus and 
even Buthns occitcinus , all thermophylous species typical 
from dry conditions. The lower height of Puerto de Tor- 
navacas (1010m) makes me to suspect that/, cyreni is not 
there and the record is possibly a mistake. One possibil¬ 
ity is that it was from the higher neighboring mountains. 

An account about the pattern and coloration of I. cyreni 
is in Arribas (1996) and it is especially detailed for the 
main Sistema Central massifs in Perez-Mellado et al. 
(1993). Both colour as well as the dark pattern, seem to 
be selected in accordance to the substrate characteristics. 
Overall, the background colour is brown in juveniles and 
subadult specimens, changing in different percentages to 
green in adult specimens (more frequently in males and 
becoming more vividly linked to reproductive processes). 
In populations inhabiting rocks (plenty of Rhizocarpon gr. 
geographicnm lichens) as in Gredos and the upper parts 
of Guadarrama (Penalara) green adults are more frequent 
(both males and big females). When living in rocky talus 
with sands and bare ground (as in Navacerrada area) 
brownish adult specimens are more frequent (Arribas 
1999b: Figs 9 and 10). 

Concerning the reticulate pattern, it also varies in a dif¬ 
ferent degree among the different populations depending 
on the substrates inhabited. Juveniles and subadults have 
temporal uniform or reticulate bands that coalesce during 
growth with dorsal spots (more frequently in males) giv¬ 
ing reticulated-likc patterns. Usually, there is a relation¬ 
ship between the size of the granite phanerocrystals (the 
granite-rock spotting) and the habitus of the lizards liv¬ 
ing on it. Lizards living on rocks that present large crys¬ 
tals (as for instance the Bejar ones) are more reticulated 
than specimens living in places in which the rocks pres¬ 
ent smaller crystals (and thus finely spotted). Coloration 
accounts described in Perez Mellado et al. (1993) are fair¬ 
ly precise, especially for Gredos specimens. Concerning 
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Bejar (=Candelario in Perez Mcllado ct al. [ 1993]) the de¬ 
scription should be corrected as, although it is true that 
very old specimens arc fairly reticulated, especially males, 
females more frequently have two paravertebral rows of 
distinctive spots (photo 2), as in females of other popu¬ 
lations and in mid-grow specimens of both sexes in all lo¬ 
calities. The statement that “the most common background 
colour of the back and flanks is greenish or bluish” prob¬ 
ably is true for fully adult specimens during the breeding 
period, but in July and August only some big males con¬ 
serve greenish tones, appearing even fully adult females 
more or less brownish. 

Concluding: a) 7. c. cyreni is a recently differentiated but 
well diagnosable (morpho)subspecies, a case paralleling 
the relationship of 7. monticola monticola from Serra da 
Estrela (Portugal) with repect to 7. /;/. cantcibricci (from 
Galicia and Cantabrian Mts.), in this later case, with no 
genetic differences, but with singular morphological traits 
that distinguishes it from other 7. monticola in a multivari¬ 
ate analysis (Arribas & Carranza 2004; Arribas et al. 
2006). Highly variable nuclear markers as for instance in- 
trons or microsatellites may help to clarify definitively the 
status of these well diagnosable (morpho)subspecies. 

b) Despite the mtDNA differences of the Bejar specimens 
(one or two nucleotides), their morphology is largely 
equivalent to 7. c. castiliana. Lacking data from other Be¬ 
jar populations and genetic nuclear markers, I assume that 
these morphological similarities reflect the presence of a 
very recent gene flow with other neighbouring popula¬ 
tions. The Bejar populations are however outstanding by 
their female body elongation (up to near 1 cm larger), con¬ 
serving proportionately shorter limbs, which can be a strat¬ 
egy to increase clutch size. 

c) The Sierras de Avila populations are very similar; the 
closer among all the populations compared. One of them, 
the Sierra de Villafranca, is together with Bejar, the most 
connected sample, and it is, from a morphological point 
of view, at the root of the species expansion. Both also 
present the unique slightly variant haplotypes. All they 
should be considered as 7. c. castiliana by their closer iden¬ 
tity with Gredos. 
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APPENDIX I 

Neotype designation for Iberolacerta cyreni (Muller & Hellmich, 
1937) 

The type series of Iberolacerta cyreni (Lacerta monticola cyreni 
Muller & Hellmich, 1937) included 66 specimens (not only from 
Guadarrama, the species' type locality, but also some paratypes 
from Credos) formerly deposited in the Zool. Staatssamlung 
(Herpet. Samml.) Miinchen. In fact, two syntype specimens, 
male and female, respectively numbered "ZSM (SLM) 2329 a” 
and “ZSM (SLM) 2329 b”, were originally considered as types 
labelled "Guadarrama, Puerto de Navacerrada. W. Hellmich”. 
(Muller & Hellmich, 1937). 

Although it seems that all the original type series was destroyed 
during the Second World War (Franzen & Glaw, 2007), and due 
to the fact that in this type series there were included some Gre- 
dos specimens (today part of another subspecies; I. c. castiliana ), 
and also that there was early confusion about the I. monticola 
type locality which lead to the description of a new taxon as Lac¬ 
erta esfrellensis Cyren, 1928 (Arribas 2008), I design a new type 
specimen (neotype) to fix unequivocally the type locality against 
any contingence (as could be the highly improbable apparition 
of any "surviving” original Gredos paratype). 

I designate here as NEOTYPE for the species a specimen from 
the Museo Nacional de Ciencias Naturales (Madrid) (MNCN n. 
39934) (Fig. 5). 

A male labelled as follows: 

Left hindleg: (white label, Typewriter letter) MNCN (anverse), 
39934 (reverse). 

Right hindleg: (white label, pencil handwritten) “Pto. de 
Cotos-Pto de Navacerrada. Srra. de Guadarrama ( )[blank in¬ 
side parenthesis], 21-IV-84, 18,15 h. Sol. Pedriza en pinar con 
nieve. Ps=8.5 gr.” (no collector’s data). 

Left foreleg: (White label, ink handwritten) Neotypus. O. Ar¬ 
ribas designatio (anverse), “ Lacerta monticola cyreni Muller & 
Hellmich, 1937” (= Iberolacerta cyreni) (reverse). 

Right foreleg: (Red plastic label, Dymo® lettering) NEOTYPUS. 

Neotype description (Fig. 5): 

Biometry: Adult male with snout-vent length of 66.85 mm. Tail 
126 mm (intact). Forelimb length 23.34 mm. Hindlimb length 
34.52 mm. Pileus length 16.4 mm. Pileus width 8.2 mm. Pari¬ 
etal legth 5.5 mm. Masseteric widest diameter 2.71 mm. Tym¬ 
panic widest diameter 1.94 mm. Anal plate width 5.04 mm. Anal 
length 3.19 mm. FLL/SVL (relative forelimb length): 0.349. 
HLL/SVL (relative hindlimb length): 0.5163. PL/PW (pileus 


shape): 2.003. DM/PaL (relative masseteric plate size): 0.491. 
DT/PaL (relative tympanic size): 0.352. AL/AW (anal plate sur¬ 
face): 0.6329. AS/SVL (relative anal plate size in respect to to¬ 
tal legth): 5.998. 

Scalation: Number of supraciliary granules: 9 (rigth) and 11 
(left). Supralabials: 5 (both sides). Sublabials: 6 (right side) and 
7 (left side). Submaxillars: 6 (both sides). Gularia: 25. Collar- 
ia: 9. Dorsalia: 53. Ventralia: 26. Femoral pores 19 (right) and 
18 (left). Lamellae: 25. Circumanal Plates: 8. Rostral in full con¬ 
tact with frontonasal. Supranasal separated from first loreal. One 
postnasal (in both sides). First Postocular separated from Pari¬ 
etal plate. Alternate wide and narrow scale rings in the tail. Twen¬ 
ty-six scales across one of these rings. 

Coloration: (in alcohol). Dorsal tract and pileus brown (proba¬ 
bly also in life), densely spotted of medium-sized black spots 
that in the middle of the dorsum nearly form transverse bands 
and connect the two temporal bands. Temporal bands reticulat¬ 
ed fairly dark (black-brown) with traces of clear occelli (bare¬ 
ly visible) inside, also connecting with more light reticulated with 
infratemporal band (barely discernible). No blue axillar occel¬ 
li. Traces of blue points in the outermost ventral scales. Only 
the outermost ventral scale ranges are clearly black spotted. Bel¬ 
ly light bluish or white-bluish. 


APPENDIX II 

Barriers and high mountain passes (among parentheses) between 
the different /. cyreni populations. All these intermediate areas 
are at present apparently devoid of /. cyreni , thus constituting 
these OTUs discrete populations: 

BEJAR-GREDOS: no barrier (Puerto de Tornavacas, 1275m). 

BEJAR-VILLAFRANCA: Tormes River Valley (no pass). 

GREDOS-VILLAFRANCA: Tonnes River Valley (Collado de 
Cepegato, 1550m). 

GREDOS-S a AVILA[Serrota+Paramera]: Alberche River Valley 
(no passes). 

SERROTA-PARAMERA: no barrier (Puerto de Menga, 
1566m). 

VILLAFRANCA-SERROTA: no barrier (Puerto de Chia, 
1663m). 

S a AVILA (as a whole) - GUADARRAMA: no barrier (Puerto 
del Boqueron, 1315m). 

GREDOS-GUADARRAMA: Alberche River Valley (no pass) 
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Table 2. As in table 1, but for/, cyreni females. 
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Table 3. Analysis of Similarity (ANOSIM) results (with 1000 randomizations). Males above diagonal and females below. The 
number is the test probability results among each two concerned populations (significant results in bold). 


FEM \ MAL 

BEJAR 

GREDOS 

GUAD. 

MIJARES 

SERR-PAR 

VILLAFR. 

BEJAR 

— 

0.058 

0.001 

0.644 

0.138 

0.064 

GREDOS 

0.021 

- 

0.001 

0.691 

0.867 

0.198 

GUAD. 

0.001 

0.001 

- 

0.097 

0.128 

0.001 

MI J ARES 

0.53 

0.039 

0.327 

- 

0.224 

0.513 

SERR-PAR 

0.127 

0.025 

0.001 

0.217 

- 

0.86 
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